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Abstract. We consider an initial- and Dirichlet boundary- value problem for a linear Cahn-Hilliard- 
Cook equation, in one space dimension, forced by the space derivative of a space-time white noise. First, 
we propose an approximate regularized stochastic parabolic problem discretizing the noise using linear 
splines. Then fully-discrete approximations to the solution of the regularized problem are constructed 
using, for the discretization in space, a Galerkin finite element method based on H 2 — piecewise polyno- 
mials, and, for time-stepping, the Backward Euler method. Finally, we derive strong a priori estimates 
for the modeling error and for the numerical approximation error to the solution of the regularized 
problem. 



1. Introduction 

Let T > 0, D = (0, 1) and (0,-7-", P) be a complete probability space. Then we consider the following 
model initial- and Dirichlet boundary- value problem for a linear Cahn-Hilliard-Cook equation: find a 
stochastic function u : [0, T] xD->R such that 

d t u + d 4 x u + fi d 2 x u = d x W(t, x) V (t, x) G (0, T] x D, 

(1.1) d 2 x m u(t,-)\ gD =0 V* 6(0,11, m = 0,l, 

u(0,x)=0 Vie£», 

a.s. in fi, where W denotes a space-time white noise on [0, T] x D (see, e.g., [23], [H]) and /x is a real 
constant for which there exists n £ N such that 

(1.2) (k - l) 2 vr 2 < n < k 2 tt 2 , 

where N is the set of all positive integers. The above stochastic partial differential equation combines two 
independent characteristics. On the one hand it corresponds to the linearization of the Cahn-Hilliard- 
Cook equation around a homogeneous initial state, in the spinodal region, that governs the dynamics 
of spinodal decomposition in metal alloys; see e.g. [4], and references therein. On the other hand the 
forcing noise is a derivative of a space-time white noise that physically arises in generalized Cahn-Hilliard 
equations, which are equations of conservative type describing the evolution of an order parameter in 
phase transitions (see [TO]: cf. [12]. |2). [19]). 

The mild solution of the problem above (cf. [B]) is given by the formula 

(1.3) u(t,x)= f f 9(t-8;x,y)dW(a,y), 

JO J D 

where 

oo 

(1.4) ^(t;x } y) = -Y / e- xl{xl -^ t e k (x)e / k (y) V (t, x, y) € (0, T] x D x D, 

fe=i 
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with Afe := kir for k 6 N, and £fc(z) := ypL sin(Afe z) for z 6 D and t e N. Observe that a;, y) = 
-d y G(t;x,y), where G(t;x,y) = Y%Li e ~ A * (A2fe_Al)t £fc(a;) e fe (y) for all (i,a:,y) G (Q,T] x 75 x U, is the 
space-time Green kernel of the corresponding deterministic parabolic problem: find a deterministic func- 
tion u:[0,r]xl)4l such that 

dtW + d£w + nd%w=Q V (t, x) G (0, T] x D, 

(1.5) 3?MV)L = V te(0,n m = 0,l, 

ui(0,x) = wo(x) Viei). 

The goal of the paper at hand is to propose and analyze a methodology of constructing finite element 
approximations to u. 

1.1. The regularized problem. Our first step is to construct below an approximate to (jl.ip regularized 
problem getting inspiration from the work [1] for the stochastic heat equation with additive space-time 
white noise (cf. [14]. [15]). 

Let N+ G N, At := J+ e N and Aa; :— j-. Then, consider a partition of the interval [0,T] 
with nodes (i n ) n f_o and a partition of D with nodes (xj) J ^ , given by t n := n At for n = 0, . . . , N± and 
Xj := j Ax for j = 0, J*. Also, set T n := (t n -i,t n ) for n = l,...,iV*, and £)j := (xj_i, :Cj) for 
j = l,...,J*. _ 

First, we let 5* be the space of functions which are continuous on D and piecewise linear over the 
above specified partition of D, i.e., 

S*:={sGC(:D;R): s\ d €¥ 1 (D j ) for j = 1, . . . , J* } C H\D). 

It is well-known that dim(<S*) = J* + 1 and that the functions (i/'i)^ 1 C 5* defined by: 

:= 35 ( x i ~ : = 35 ( x ~~ a; ^-i) + ' 

^i(x) := [ - Zi-a) ^x^,*^] + (as - a*) «*(» < _ 1> «, j ] ] > i = 2, J*, 

consist the well-known hat functions basis of <S*, where, for any 4 C 1, by X A we denote the index 
function of A. Next, consider the fourth-order linear stochastic parabolic problem: 

d t u + d%u, + vdlu = d x W in (0, T] x Z), 

(1.6) dl m u{t,-)\ gD = Q V* 6(0,11, 171 = 0,1, 

u(0,x)=0 Viefl, 

a.e. in f2, where: 



ra=l 



£=1 \m=l 



, V (f , x) 6 [0, T] x D, 



G is a real, (J* + 1) x (J + + 1), symmetric and positive definite matrix with 

G l .j := (ipj,ipi) 0iD , i, j = 1, J* + l, 

and 

R n .i-= ipi(x) dW(t,x), i = l,..., J* + l, ri= 1,...,JV*. 

J T n J D 

The solution of the problem (|1.6[) , has the integral representation (see, e.g., |17j ) 



u(x,t) = j j G(t- s;x,y)d y W(s,y)dsdy 

(1.7) 



*(< - s; a;, y) VK(s, y) dsdy, V (t, a;) 6 [0, T] x TJ. 

0./D 



Remark 1.1. A simple computation verifies that G is a tridiagonal matrix with Gx,x = G Jt+W(+1 = 4p, 
Gi,i — ^-j^ for i — 2, . . . , J*, and Gi,i+x = Ah /or i = 1, . . . , J*. Since G is symmetric we have in 
addition that Gi-x.i = /or i = 2, . . . , J* + 1. 

Remark 1.2. Let T = {(n, i) : n = 1, . . . , iV*, i = 1, . . . , J* + 1}. Using the properties of the stochastic 
integral (see, e.g., [23 \), we conclude that R n ^ ~ A/"(0, At Gj,j) /or a/Z (n, j) € I. ^Mso, we observe that 
E[R n _i R n ' ,j] = /or (n,i), (n',j) 6 I wi£/i n 7^ n', and hence they are independent since they are 
Gaussian. In addition, we have that E[i? n ^ R n ,j] = AtGij for (n,i), (n,j) G 1. Thus, for a given n the 
random variables {Rn^Yi^X 1 o,re Gaussian and correlated, with correlation matrix AtG. 

1.2. The numerical method. Our second step is to construct finite element approximations of the 
solution u to the regularized problem. 

Let M G N, At :— -p., r m :— mAr for m = 0, . . . , M, and A m :— (r m _i,r m ) for m — 1, . . . ,M. Also, 
let r G {2,3}, and M£ c H 2 (D) n Hq(D) be a finite element space consisting of functions which are 
piecewise polynomials of degree at most r over a partition of -D in intervals with maximum mesh-length 
h. Then, computable fully-discrete approximations of u are constructed by using the Backward Euler 
finite element method, which first sets 

(1.8) U° h := 

and then, for m = 1, . . . , M, finds U™ G M£ such that 



(1.9) (UF-U™-\x)o.d+*T ((U?yW')o.n+ti((U?r,x)o,D = (d x W, X )o,odT 

for all x € where (•, -) D is the usual L 2 (D)— inner product. 

1.3. An overview of the paper and related references. Our analysis first focus on the estimation 
of the modeling error, i.e. the difference u — u, in terms of the discretization parameters At and Aa;. 
Indeed, working with the integral representation of u and u, we obtain (see Theorem 13. 1[) 
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(1.10) max I J f J \u(t,x) -u(t,x)\ 2 dx^j dP^ 2 < C mc (e~i Ax i *- t + Ati), Vee(0,|], 

where G mc is a positive constant that is independent of Ax, At and e. Next target in our analysis, is to 
provide the fully discrete approximations of u defined in Section ll.2l with a convergence result, which is 
achieved by proving the following strong error estimate (see Theorem 15. 3|) 



(1.11) ^m^^J (J \Uj^(x) -u(r m ,x)\ 2 dx^j dP^ 2 < G nc ( e l 1 Ar^ £l + e 2 ^"W-ft ) , 

for all €1 G (0, |] and e 2 G (0, v(r)] with v(2) — i and ^(3) = \, where G ne is a positive constant 
independent of ex, e 2 , At, h, Ax and At. To get the error estimate (jl.lip we use as an auxilliary tool 
the Backward-Euler time-discrete approximations of u which are defined in Section [H Thus, we can see 
the numerical approximation error as a sum of two types of error: the time-discretization error and the 
space- discretization error. The time- discretization error is the approximation error of the Backward Euler 
time-discrete approximations which is estimated in Theorem 14.21 while the space- discretization error is 
the error of approximating the Backward Euler time-discrete approximations by the Backward Euler 
finite element approximations, which is estimated in Proposition 15.21 

Let us expose some related bibliography. The work [18] contains a general convergence analysis for a 
class of time-discrete approximations to the solution of stochastic parabolic problems, the assumptions of 
which may cover problem However, the approach we adopt here is different since first we introduce 

a space-time discretization of the noise and then we analyze time-discrete approximations to the solution. 
We would like to note that we are not aware of another work providing a rigorous convergence analysis 
for fully discrete finite element approximations to a stochastic parabolic equation forced by the space 
derivative of a space-time white noise. We refer the reader to our previous work [TJ], [TS] and to [TB] for 
the construction and the convergence analysis of Backward Euler finite element approximations of the 
solution to the problem (jl.lj) when fi = and an additive space-time white noise W is forced instead of 
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d x W . Finally, we refer the reader to [8], [T], [13], [3J, [22] and [23] for the analysis of the finite element 
method for second order stochastic parabolic problems forced by an additive space-time white noise. 

We close the section by an overview of the paper. Section [5] introduces notation, and recalls or proves 
several results often used in the paper. Section |3j is dedicated to the estimation of the modeling error. 
Section [4] defines the Backward Euler time-discrete approximations of u and analyzes its convergence. 
Section [5] contains the error analysis for the Backward Euler fully-discrete approximations of u. 

2. Notation and Preliminaries 

2.1. Function spaces and operators. Let / C I be a bounded interval. We denote by L 2 (I) the 
space of the Lebesgue measurable functions which are square integrable on / with respect to Lebesgue's 

measure dx, provided with the standard norm ||ff|| 1 '■= (I \g{x)\ 2 dx) 2 for g £ L 2 (I). The standard inner 

product in L 2 (I) that produces the norm || • || J is written as (-,-)o,ij i- e -j {9ii92)o,i := /, 9i{ x )92{x) dx 

for gi, <?2 £ L 2 (I). Let No be the set of the nonnegative integers. For s £ No, H S (I) will be the Sobolev 

space of functions having generalized derivatives up to order s in the space L 2 (I), and by |j • \\ SiI its usual 

i 

norm, i.e. := (J2e=o ll^ £ 5llo,i) 2 f° r 9 ^ S {I). Also, by Hq(I) we denote the subspace of H X {I) 

consisting of functions which vanish at the endpoints of I in the sense of trace. We note that in Hq(I) 
the, well-known, Poincare-Friedrich inequality holds, i.e., there exists a nonegative constant C PF such 
that 

(2.1) IMk, < CW||00|| o ,, VgeH^I). 

The sequence of pairs ( (X 2 ,, £&) is a solution to the eigenvalue/eigenfunction problem: find 
nonzero <p G H 2 (D) n Hq{D) and a £ R such that ~d 2 ip = a ip in D. Since (£k)kLi is a complete 
(•, -) iD — orthonormal system in L 2 (D), for s £ R, a subspace V S (D) of L 2 (D) is defined by 

V S (D) := L £ L 2 (D) : £ X 2s (v,s k )l D < oo j 

i 

which is provided with the norm \\v\\ v s := ( J2kLi ^l? ( w > e k)1 D ) 2 ^ ' v e V S (D). For s > 0, the pair 
(V S {D),\\ ■ \\ v s) is a complete subspace of L 2 {D) and we set (H S (L>), || • := (V s (£>),|| • || v ,). For 
s < 0, we define (H S (D), || • ||h=) as the completion of (V s (D), || ■ || v s), or, equivalently, as the dual of 
(H~ S (L>), || ■ \\ A -s). Let m G N . It is well-known (see [21]) that 

(2.2) H m (L>) = { v £ H m (D) : 3 2l v \ BD = if < i < f } 
and there exist positive constants C m . A and C m . B such that 

(2.3) C m , A \\v\\ m , D < ||v|| A m < C m , B ||7j|| m , D , V« G H m (yD). 
Also, we define on L 2 (D) the negative norm || • ||_ m , D by 

\H-m, D :=s U p{^f^: pe^fl) and (p^o}, Vi> G L 2 (D), 

for which, using (|2.3p . it is easy to conclude that there exists a constant C_ m > such that 

(2.4) ||»||- m . fl < C_ ro ||«|U-», Vv£L 2 {D). 

Let L2 = (L 2 (D), (•, -) ,r>) and £(L2) be the space of linear, bounded operators from L2 to L2. We 

say that, an operator T G £(L 2 ) is Hilbert-Schmidt, when ||r|| H s : = (S^i l|r(efe)||o D ) 2 < +00, where 
||r|| H s is the so called Hilbert-Schmidt norm of T. We note that the quantity ||r|| H s does not change 
when we replace (e/ s )^_ 1 by another complete orthonormal system of L2, as it is the sequence (<fk)kLo 
with ipo(z) := 1 and tpk{x) := \/2 cos(Afe z) for k £ N and z £ D. It is well known (see, e.g., [7]) that an 
operator T £ £(1/2) is Hilbert-Schmidt iff there exists a measurable function g : D x Z? — > R such that 
(r(u))(-) = J D g(-,y)v(y)dy for u G L 2 (D), and then, it holds that 

(2.5) ||r|| HS = ^ J g 2 (x,y)dxdy 
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Let £ H s(L2) be the set of Hilbert Schmidt operators of £(L 2 ) and $ : [0, T] — > £ H s(L2)- Also, for 
a random variable X, let ELY] be its expected value, i.e., K[X] := j XdP, Then, the Ito isometry 
property for stochastic integrals, which we will use often in the paper, reads 



(2.6) 









2 




E 




/ $ dW 




-L 






Jo 


O.D 





l$(t)lll s dt. 



Let II : L 2 ((0,T) x5)-> L 2 ((0,T) x D) be a projection operator defined by 



g(s,y)ipi(y) dsdy tpi(x), V(t,x) G T n x D, 



(2.7) fij(M):=ii J f / / 

for n = 1, . . . , N* and for g G L 2 ((0, T) x L>), for which holds that 



(ihj) 2 dx* 



Jd 



< 



g A dxdt) , V 5 G Z/((0,T) x £>). 



Now, in the lemma below, we relate the stochastic integral of the projection LI of a deterministic func- 
tion to its space-time L 2 — inner product with the discrete space-time white noise kernel W defined in 
Section 0(cf. Lemma 2.1 in [H]). 



Lemma 2.1. For g G L 2 ((0,T) x D), it holds that 



(2.9) 



Ilg(t,x)dW(t,x) = 



J D 



W(s,y)g(s,y) dsdy. 



J D 



Proof. To obtain (|2.9[) we work, using ()2.7j) and the properties of the stochastic integral, as follows: 
/ / Ug(t, x) dW(t, x)=i- t J2J2J2 G ll / 9(s, y) Mv) dsdy R na 

JO Jd n=l i=1 £ = 1 \JT n xD J 

=S E / 9(s, y) (f £ G#lMl/) J^i ] 



71=1 



_-, ^0 •/£> 



1=1 £=1 



fif(s,y) W(s,y) dsdy. 



We close this section by observing that: if c* > 0, then 



□ 



(2.10) 



fe=i 



-(l+c„e) 



< 



(^)i Ve G (0,2], 



and if (H, (•, -) w ) is a real inner product space, then 

(2.H) (5 - > \ i(g,g)n- {v,v)h}, Vg,veH. 

2.2. Linear elliptic and parabolic operators. Let us define the elliptic differential operators A B , A B : 
H 4 (D) — > L 2 (D) by A B v := d A v + [id 2 v and A B v :— A B v + fi 2 v for v G H 4 (D), and consider the 
corresponding Dirichlet fourth-order two-point boundary value problems: given / G L 2 (D) find v B , 
v B G H 4 (L>) such that 

(2.12) A B v B = / in D 
and 

(2.13) A B v B = f in D. 

5 



Assumption (|1.2[) yields that when k — 1 or k > 2 and fx ^ A^_ l5 the operator A B is invertible and thus 
the problem (|2.12l) is well-posed. However, the problem (|2.13l) is always well-posed. Letting T B , T B : 
L 2 (D) -t H 4 (L>) be the solution operator of ([2~T2l and ([2^15]) . respectively, i.e. T B f := A" 1 / = v B and 
T B f := A^ 1 / = v B , it is easy to verify that 

oo oo 

(2.14) T B f = Y,M^ £k and f -/ = E#fe? e *. V/GL 2 P), 

fc=l k-l 

and 

(2.15) ||T B /|| ra>D + ||f B /|| m , D < C R , m ||/|| m _ 4>Dj V/ g fpn«{o,™-4}( jD ) i vm G N 0> 

where C R m is a positive constant which is independent of / but depends on the D and m. Observing 
that 

(T B v 1 ,v 2 ) ,n = (vi,f B v 2 ) ,o, Vv!,v 2 G L 2 (D), 
and in view ([27T4j) . the map 7 B : £ 2 (£>) x i 2 (L>) -> K defined by 

j B (v,w) — (T B v,w) 0}D Vv, u> G L 2 (D), 

is an inner product on L 2 (D). 

Let (<S(t)iuo)te[o,T] be the standard semigroup notation for the solution u> of (| 1 . 5[) . Then, the following 
a priori bounds hold (see Appendix fA")) : for f G No, /? > and p > 0, there exists a constant Cp t p ^ > 
such that: 

(2-16) J^iT-tafWdfSi^woW^dT^ Cp A ^ T \\w \\l P+Ae ^ 2 

forall w G HP+ 4£ - 2 ^ 2 (L>) and t a , t b G [0,T] with i 6 > t a . 

2.3. Discrete spaces and operators. For r G {2,3}, let M£ C Hq(D) PiH 2 (D) be a finite element 
space consisting of functions which are piecewise polynomials of degree at most r over a partition of D in 
intervals with maximum mesh- length h. It is well-known (cf., e.g., [5]) that the following approximation 
property holds: 

(2.17) inf ||«-x|| 9 , D < CW^-'HUi,^ Vve H s+1 (D)nH^(D), Vse{2,r}, 

where C FM ^ r is a positive constant that depends on r and is independent of h and v. Then, we define the 
discrete elliptic operators A B _ h , A Bih : M£ — > M£ by 

(2-18) (A B ^,x)o,o:=(d 2 ip,d 2 x)o,o+fi(d 2 <p,x)a,n, V^, X €M£, 

and 

(2.19) A Bth <p := A Bih ip + n 2 <p, Vy>GM£. 

Also, let Ph '■ L 2 (D) — > M£ be the usual L 2 (D)— projection operator onto M£ for which it holds that 

(Phf,x)o,» = (f,x)o,», VX6MJ, \ffeL 2 (D). 
A finite element approximation S BJ , G M£ of the solution v B of (|2.13|) is defined by the requirement 

(2.20) A B , h v B>h = P h f, 
where the operator A B h is invertible since 

(2-21) (A B , hX , X)o.o > \ ( ||9 2 X|| 2 , D + M 2 llxll 2 D ) , V X G M£. 

Thus, we denote by T Bfl : L 2 {D) -> M£ the solution operator of (|2.20p . i.e. 

T B , h f:=v Bih =A-\P h f, VfeL 2 (D). 
Next, we derive an L 2 (D) error estimate for the finite element method (|2.20j) . 



Proposition 2.1. Let r g {2,3}. Then we have 

(h 4 \\f\\ , Dl r = 3, 

(2.22) ||T B / - f B<h /|| , fl <clh 3 ||/||_ liD , r = 3, V/ e £ 2 (£), 

U 2 r = 2, 

where C is a positive constant independent of h and f. 

Proof. Let / e L 2 (D), e = f B f - f B h f and v = f B e. To simplify the notation we define B : H 2 (D) x 
H 2 (D) — > R by B(v, w) := (d 2 v, d 2 w) D + fi (d 2 v, w) , D + fi 2 (v, w) , D for v, w e H 2 (D). It is easily seen 
that 

(2 23) B(v,w)< V^l + M) (\\d 2 v\\l»+» 2 \Ml D y \Hk» Vv lW eH 2 (D), 

B(v,v)> k[\\dMln+V 2 \\v\\ln] VveH 2 (D). 

Later in the proof we shall use the symbol C for a generic constant that is independent of h and /, and 
may changes value from one line to the other. 

First, we observe that ||e|| 2 D = B(e,v). Then, we use the Galerkin orthogonality to get 

\\e\\l D =B(e,v- X ), V X eA4 r , 

which, along with (|2.23|) . leads to 

(2-24) || e || 2 D < C (\\d 2 e\\l D +^\\e\\l D y ini l\\v - x\U.». 

Using again (|2.23p and the Galerkin orthogonality, we obtain 
\\d 2 e\\l D +^\\e\\l D <2B(e,e) 

<2B(e,f B f- X ) 

<C (\\d 2 e\\l D+f i 2 \\e\\ 2 0D y \\T B f- X \U,n, V^M[, 

which yields that 

(2-25) {\\d 2 e\\l D +^ 2 \\e\\l D y < C inf ||T B / - X || 3 , D . 

Combining (j2~24l) . ([2~25) and (j2~T71) , we arrive at 

||e|| 2 , D <C inf ' \\T B f- X \Un inf \\v- X \k» 
(2.26) xeMh x ^ Mh 

< C h s+s '- 2 ||T B /|| s+ltD \\T B e\\ s , +1>D , V«, s' € {2, r}. 

Let r = 2. We use (j2~26l) and (|2~T5)) to get 

IMI 2 , D <c^ 2 ||T B /|| 3 , D ||T Be || 3 , D 

<Ch 2 \\f\U D \\e\U D 
<Ch 2 \\f\\^ D \\e\\ 0:D , 

from which we conclude (|2.22[) for r = 2. 

Let r = 3. We use (|2~26l) with s' = 3 and (|2~T5j) to obtain 

||e|| 2 , D <C/i s+1 ||T B /|| s+1<D ||T B e|| 4 , D 

<Ch s+1 \\f\\ s _ 3 , D ||e|| D>D , s = 2,3, 

from which we conclude (|2.22|) for r = 3. □ 



Let 7 B , h : L 2 (D) x L 2 (D) -> R be defined by 

7a. h (/,0) - (T B ,J,g) 0lD V/,.9 e L 2 (£). 
Then, as a simple consequence of (|2.2ip . the following inequality holds 

(2.27) %MJ)> ^ (pHT B ,J)f , D +» 2 \\T B ,J\\l D ), VfeL 2 (D). 

Thus, observing that 

(T Bi J,g) 0>D = (f,T B , h g) , D , Vf,g € L 2 (D), 
and using (|2.27p . we easily conclude that j B h is an inner product in L 2 (D). We close this section with 
the following useful lemma. 

Lemma 2.2. There exists a positive constant C > such that 

(2-28) %MJ)< C\\ff_ 2tD , VfeL 2 (D). 

Proof. Let / e L 2 {D), ip = f B f and ip h = f BJ J. Then, we have 

(T B ,hf,f)o.n = (A B V>, 1ph)o,n 

(2.29) = (9V, 9Vft)o, D + A* ^)o, D + W>, Ww 

< \ ( II#V|| 2 , D + M 2 ll^llo 2 d ) + e ( ll^lliU + M 2 \\M\lo ), Ve > 0. 
Setting e = i in (|2.29p and then combining it with (|2.27p . we obtain 

(2-30) ll^l| 2 , D +M 2 H^llo, D < M{\\d 2 iP\\l D +V 2 Mlo)- 

Finally, (j2~29| with e = |, (j230l) and (l2~T5)) yield 

7b,,(/,/)<8 (ii^llL + ^IHliU) 

<8(l + M 2 )||f B /|| 2 D 
<8(1 + ^ 2 )C^ 2 ||/|| 2 2 , D . 

Thus, we arrived at (|2~28)) . □ 

3. An Estimate for the Modeling Error 

In this section, we estimate the modeling error in terms of At and Ax (cf. Theorem 3.1 in [14]). 

Theorem 3.1. Let u be the solution of and u be the solution of (|1.6p . Then, there exists a real 

constant C > 0, independent of At and Ax, such that 



(3.1) 



max(E [||u-u|| 2 D ] ) 2 < C |u (At) At * + e"= As*" 



Vee(0,i], 



w/iere w Q (Af) := V 1 + At?. 

Proof. Using (|1.3I) . (|1.7p and Lemma |2~TI we conclude that 



(3.2) u(t,x)-u(t,x) = [X (0}t) (s)*(t-s;x,y)-*(t,x;s,y)]dW(s,y), V (t, x) e [0, T] x D, 



ds', V(s,y)eT„x 



where f :(0,T)xD-) L 2 ((0, T) x £>) is given by 



for n = 1 , . . . , iV* . 



Let 6 := {E - u\\ 2 D ] } 2 and i e (0, T]. Using (O and Ito isometry we obtain 



e(t) 



#(o,t) («) *(* - a; x, y) - x; s, y) 



dxdyds 



Now, we introduce the splitting 

(3.3) 

where 

e.(*) :={e 



.71=1 



_ 1 J DJ D J Ti 



e(t) < e A (t) + e fl (t), 



dxdyds 



and 



e*(*) := E 



,rt=l 



D J D J T n 



X m (s)^(t-s;x,y)-^ / # (0)t) (s') ¥(i - a';a:,i/) da' 



dxdyds 



Also, to simplify the notation in the rest of the proof, we set \i k := \\ (A 2 . — fj,) for k e N, and use the 
symbol C to denote a generic constant that is independent of At and Ax and may changes value from 
one line to the other. 

• Estimation of Q A (t); Using (|1.4|) and the (•, -) D — orthogonality of {£k)kLii we nave 

-i 2 



<£(*) = &E 



#(„,,)(*') *(t-y ; x,l/)- E G-, 1 (*(t- S ';x,-),^(-)) , D ^(2/) 



ds' 



dydx 



At E 



El / ^(^Oe^^ds' / - E G *~i (4>^)o, D &(v)) rfy 
.fc=i ' Jd 

from which, using the Cauchy-Schwarz inequality, follows that 

K OO 

(3.4) e 2 A (t)<Y / Mt)B k + E Mt)B k , 

k=l fc=«+l 

where 

Afc(i) :=2A^ f e- 2 ^^~ s ">ds', 
Jo 

2 



B,, 



k ■= / (vfe(y) - E (fk^i) ,D ip%{yij dy. 



First, we observe that 



(3.5) 



Bit < max sup \<Pk(x) - <Pk(y)\ 



< min{l, AfcAa;} 

< min{l,(V2A fc A2;) e }, V6»e[0,l], VfceN. 



A fe (i)<i^^ 



Next, we use (|1.2[) . to obtain 
(3.6) 

Thus, from p.4[) , (|3.5p and (|3.6p . we conclude that 



< I iw^ vfc>« + i. 



el(t)< C ( (Ax) 2 E^ + (Aa;) 2e E ^ 

k=K+l h 



k=l 



which yields 
(3.7) 



e A (t)<c(Axf [£; v^[o,i; 



_ i + 2(i _e) 
.fe=l A fc 
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Estimation of 9 B (t): For t G (0,T], let N(t) := min{ I G N : 1 < I < N* and t < t £ } and 

' T„, if n< #(t) 



r n (t) :=T B n(o,t) 



n = l,...,iV(t). 



(t^w-nt), if n = JV(t) 
Thus, using (|1.4|) and the (•, -) 0)D — orthogonality of {ek)kLi and {^kj^Li as follows 



e s(*) - (SF X! 

n=l 

= wE 

71=1 

we conclude that 



(3.8) 
where 



D J D J T n 



x (o,t)(s) ®(t - s; x, y) ~ Afot) (s) ®(t - s'; x, y) 



ds 



dxdyds 



DJ D J T n 



J2^e k (x)My) / k |i) ( S )e-"»( t - s )-% ) ( S ')e-'"('- ! ') 



k=l 



ds' 



dxdyds 



N(t) 



© 2 B w < E A M w E *«w > 



(Ar (0lt )(a) e^(*- s ) - A- (0 , t) ( S ') e^'^) ds' 



ds. 



Let k £ N and n G {1, . . . , N(t) - 1}. Then, we have 

Mfce -/*(t-r) drds'Y ds 

max{s' ,s} 



Kit) - ( 



T n J T n J S 



< 



II k 



e -^(t-r) dTds >) dg 



< 2 



(/X, 



•) a *, +2 / ( 



H k e-^ k{t ~ T) drds') ds + 2 

<2£t(f f Hke'^-^ drds'Y + 2{At) 2 f (j fi k e '^ k{t ^ T) dr) 2 ds 
from which, after using the Cauchy-Schwarz inequality, we arrive at 



ds. 



(3.9) **(*)< 4(At) 2 ^ (jf we-^Njr) 
For k < k, we use (|3.9j) to get 

(3.10) **(f) < 4 max ( M ) 2 (At) 5 . 

i<fc<^ 

For k > k + 1 , we use (|3.9|) to have 

**(«)< 4 (At) 2 J ( e -»(*- S ) - e -»(*-*n-l)^ dfi 

(3.11) < 4(At) 2 (l-e-^ At ) 2 / 

< 2 (At) 2 (l-e - "* 



Summing with respect to n, and using (|3.9[) . p.lO|) and p. lip , we obtain 



JV(t)-l 



(3.12) 



(At) 



^ E <(*) - c { (i- 



n=l 



k < k, 

, ft > AC + 1 
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Considering, now, the case n = N(t), we have 



(3.13) 
with 

For fc < k, we obtain 
(3.14) 

For A; > k + 1 , we have 



1 *| w (t)< CAt. 



(At) 7 ^iv(t)^ 



1 - e 



< AW j_ _ -/.* At *2 



and 



"I 2 



jv(*)-i L ''JVC*)-! 



in- 



. e-v k{f ~ T) drds' + At e ^ fc(t " s) 



ds 



< 2 



< 2 



N(t) — 1 L JV(*)-1 



(is 



N(t)-1 



t (■max{s,s'} 



(At) 2 



(is 



1 - e 



— ( t— ty 



JV(t)— l ^(t)-! 



(At) 2 



( 1 - e ~ 2Alfc At ) 



N(t)-1 
t 



< 8 (At) 2 

< 8 (At) 2 

"-jv(t)-l 

which, along with (|3.13p , gives 



. e -M*(*-r) dr 



«(t)-i 



< 5^- ( 1 - e 



D -2^fc At \ i At 



1 - e 



-PfcAt ' 



«(t) - " Mfc v , ■ fj,- 

Since the mean value theorem yields: 1 — e ~ fJ - kAt < /i^ At, the above inequality takes the form 



(3.15) 



1 ,Tffe 



< 6 



l_ e -2(*fcAt 



(AFp **W - - Mfc 

Combining ([3T5]l , ([3~T2]) . (j3~14l and (l3~T5)l we obtain 



(3.16) 



9 2 (t)<C 



<C 



At+ £ A 2 fe I 

oo 



-on Af At 



fe=l 
11 



with c = 2 ^+iyi ■ To get a convergence estimate we have to exploit the way the series depends on At 
in the above relation: 



k=l 



< 



l-e~°o 1 



l- e - c o ■■ 



■ dx 



(3.17) 



< C 

< C 

< C 



1 



-CQ 7T 4 At ' 



At 



x 2 e -c X* 7T At 



At + (At) 
(At) 3 + 1 



y 2 e~ 2y dy 



(Aty 



Using the bounds (|3.16p and p,17[) we conclude that 
(3.18) 9 B (t) < C [(At)* + 1 



At*. 



The error bound (|3.1j) follows by observing that 0(0) = and combining the bounds (|3.3[) . (|3.7|l . 
dnigD and ([530]) . □ 

4. Time-Discrete Approximations 

The Backward Euler time-stepping method for problem (jl.6p specifies an approximation ?7 m of u(r m , •) 
starting by setting 

(4.1) U°:=0, 

and then, for m = 1, . . . , M, by finding U m € H 4 (Z?) such that 



(4.2) 



U m - U % 



At A B U r ' 



d x W ds a.s.. 



The method is well-defined when the differential operator Q b , At :— I + At A B : H 4 (Z?) — > L 2 (D) 
is invertible. It is easily seen that Q Bi&T is invertible when 1 + At X 2 (A 2 , — fi) ^ for fc 6 N, or 
equivalently when: k = 1 or /c > 2 and At max A 2 . (^ — A 2 .) 7^ 1. If K > 2, then it is easily seen 

that max A 2 (/i — A 2 .) < so the condition At ^- < 1 is a sufficient condition for the invertibility of 

X fc^C /*£ — X 

Qb,At- 

4.1. The Deterministic Case. The Backward Euler time-discrete approximations of the solution w to 
the deterministic problem (|1.5|) are defined as follows: first we set 

(4.3) W° := wq, 

and then, for m = 1, . . . , M, we find W m € H 4 (D) such that 



(4.4) 



W 



m — 1 



A T k R W m = 0. 



Obviously, the Backward Euler time-discrete approximations are well-defined when Q b ,/\t is invertible. 
Our next step, is to derive an error estimate in a discrete in time L\ (L 2 .) norm, taking into account that, 
in constrast to the case y, = considered in |14j . the operator is not always invertible. 

Proposition 4.1. Let (lU m )^ = o ^ e ^ e Backward Euler time-discrete approximations of the solution w 
of the problem (jl.5|) defined in (|4.3[) - (|4.4p . Also, we assume that k = 1, or n>2 and At /i 2 < j. Then, 
there exists a constant C > 0, independent of At, such that 



(4.5) 



W(T m , - J 



< C(AT) e II w ViB eH 2 (fl), V6»e[0,l]. 
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Proof. The estimate (|4.5|) will be established by interpolation, after proving it for 8 = 1 and 6 = 0. 

Let wo G H. 2 (D). According to the discussion in the begining of this section, when k = 1 or n > 2 and 
At n 2 < j, the existence and uniqueness of the time-discrete approximations (W m )m=o is secured. We 
omit the case k = 1 since then the operator A B is invertible and the proof of (|4.5[) follows moving along 
the lines of the proof of Proposition 4.1 in [T3], or alternatively moving along the lines of the proof below 
using the operator T B instead of T B . Here, we will proceed with the proof of (|4.5[) under the assumption 
At /z 2 < i , without using somewhere a possible invertibilty of A B . In the sequel, we will use the symbol 
C to denote a generic constant that is independent of At and may changes value from one line to the 
other. 

Let E m {-) := w(r m , •) - V m {-) for m = 0,.. .,M and a m := J &m [w(r m , •) - w(t, •)] dr for m = 
1, . . . , M. Then, combining (|1.5|) and (|4.4[) . we conclude that 



(4.6) T B (£ m - £'"-') + At E m = At ^ T B E m + [a m - /x T B er m J , m = 1, . . . , M. 

Now, take the L 2 (D)— inner product with of both sides of (|4.6[) . to obtain 

1b{E" 1 — E" 1 ^ 1 , E m ) D + At 1 1 £ m ||„ a = At [i 2 j B (E m , E m ) 

( 4 -7) 

+ (cr m - /i 2 T B cr m ,E m ) 0iD , m = 1, . . . , M. 

Using (j2TTj) . (|4~7)l and (j2"7T51) . we arrive at 

/ 7^™ ^-t^" 1 -^™- 1 )^^ 

(4.8) , „ 

+ (7AT- 1 || ( r m || 2 I3 , m = l,..., 

Since 2At/i 2 < 1, ([4T5]l yields 

%(E m ,E m ) < [^(E^^E^ + CAT-'Wa^lo}, m = l,...,M. 

Then, we apply a simple induction argument and use that E° — and 4Ar/x 2 < 1, to obtain 

m 

(4-9) fcl m 

^Ce^Ar-^IM^, m = l,...,M. 

Next, we use the Cauchy-Schwarz inequality to bound cr m as follows: 



cr 



2 

m\\ ,D 



<C(At) 3 [ \\d T w( S ,-)\\l D ds, m = l,...,M. 



(4.10) 

Thus, (|4TTU1) and yield 

(4.11) 7 B (£™, £T) < C (At) 2 f ™ \\d r w(s, -)\\ 2 D ds, m = 1, . . . , M. 

Jo 

Combining P~5)l . (|4~TTj) and (|4~TU|) . we have 

7 B (£; m ,i? m )-7 B (^ 1 ^ m " 1 ) + AT||i : ;™|lL <^(At) 2 / ||9 TW ( s ,.)|| 2 D d S 

( 4 - 12 ) >„ 

+ C(At) 3 / ||3 T «;( s ,.)|| 2 D d S 
Jo 

for m = 1, . . . , M. Summing with respect to m from 1 up to M and using the fact that E° = 0, (|4.12[) 
yields 

(4.13) %{E M ,E")+j^AT\\E™\\l D < C(At) 2 f \\d T w(s, .)|| % ds. 

m=l J( > 
13 



Finally, use (HT5)) and (gHH} (with /3 = 0, £ = 1, p = 0) to obtain 

(4-14) (? Ar l|£m|l °-") ^ CAr ll w «^ 2 ' 

which establishes (|4. 5[) for 9 = 1. 

First, we observe that (|4.4[) is written equivalently as 

T B (PF m - W m_1 ) + ArW m = At /x 2 T B W m , m = 1,...,M, 
from which, after taking the L 2 (D)— inner product with VF m , we obtain 

(4.15) %{W m - W m ~\ W m ) , D + At \\W m \\ 2 >D = Ar M 2 j B (W m W m ), m=l,...,M. 
Then, we combine (|2.11l) and (|4. 15)) to have 

(4.16) {l-2AT^ 2 )^ B (W m ,W m ) + 2AT\\W m \\l D < %(W m -\ IF" 1 " 1 ), m=l,...,M. 
Since 4^ 2 At < 1, (|4~TrJ)) yields that 

j B (W m ,W m ) < 1 _ 2 ^ AT 7 B (^ m - 1 ,^ m ^ 1 ) 

<e 4 ^ AT ^ B (W m -\W m - 1 ), m = 1,...,M, 
from which, applying a simple induction argument, we conclude that 

(4.17) max %(W m ,W m ) < C^ B (w ,w ). 

0<m<M 

Now, summing with respect to m from 1 up to M, and using (|4.17j) . f)4.16p yields 

M 

V AT\\W m \\ 2 D <C(f B wo,w ) , D 

(4.18) ^ 

< l|wo||-2, D IITb^oIU.d- 

Thus, using (|4T8|) . (f2~l"5j) and ([231), we obtain 



(4.19) 



In addition we have 



^2^r\\W m \\l D ) <C|K||- 2 

m=l / 



^ AT||u;(T m ,-)|| 2 D < £ / ([ d T [(T-T m ^)w 2 (T 7 x)]dT) dx 

m=l m=l"' £ ' ^A™ ' 

<£/(/ [^..) +a ( T -r^,*(r,.,-(r,.,]*)* 

m=1 JD ViA ra / 

A/ - 

< £ / [2|mT,.)||^ + (T-T m _ 1 ) 2 H^(r,-)|lUrfr 



m— 1 

< 2 



j[ [ll^-)ll^+r 2 ||^(T,.)|| 2 . D 



G?T, 



which, along with ff2~T6|) (taking (P,£,p) = (0,0,0) and (,S,£,_p) = (2, 1,0)) and jM]), yields 

1 

CM \ 2 

£ AtIHw)!! 2 , d <c|K||„- 2 . 
m=l / 

Thus, the estimate ((4~5|i for 9 = follows easily combining ([1715} and ([O0|) . □ 
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4.2. The Stochastic Case. Next theorem combines the convergence result of Proposition 14.11 with a 
discrete Duhamel's principle in order to prove a discrete in time (L 2 P (L^,)) convergence estimate for 
the time discrete approximations of u (cf. [T3], [22] ). 

Theorem 4.2. Let u be the solution of (11.61) and (C/ m )m=o ^ e ^ e time- discrete approximations defined 
by (|4.1[) - ([4.2p . Also, we assume that k = 1, or k > 2 and At/i 2 < j. Then, there exists a constant 
C > 0, independent of At, Ax and At, such that 



(4.21) 



max ( E 

1< m < m 



\U r ' 



1 < Cwi(Ar,e) At*" 6 , Ve€(0±], 



w/iere £Ji(Ar,e) := e~~- + (Ar) e (l + (Ar)i + (At)*)'. 

Proof. Let 7 : L 2 (D) -> 7 2 (7)) be the identity operator, A : 7 2 (7J) -> H 4 (7>) be the inverse el- 
liptic operator A := (7 + AtA b ) _1 which has Green function G A (x,y) — J2kLi i+a!aV-p) ' *' e ' 
A/(x) = J D G A (x,y)f(y)dy for x e D and / e 7 2 (T>). Also, we set G*(x,y) := -d y G A (x,y) = 

-Er=i i+a%-&) ' and definG * : L2(jD) by $/(2;) := L G -(^,2/)/(2/) d V for / G 7 2 (7J). 

Also, for m e N, we denote by G A *,™ the Green function of the operator A m_1 $. In the sequel, we will 
use the symbol C to denote a generic constant that is independent of At, At and Ax, and may changes 
value from one line to the other. 

Using (14.21) and a simple induction argument, we conclude that 



U m = 



J2 [ A m - j $W(T, -)dr, m = 1, ... , 



M. 



which is written, cquivalently, as follows: 



(4.22) 



U m (x) 



K. m (T;x,y)W(T,y)dydT, VieD, m = l. 



,M, 



where K, m {j;x,y) := Y^Li X ^( T ) G A^,m- j+ x(x,y), Vr€ [0,T], \/x,yeD. 

Let m G {1,...,A7} and £ m := E[||t/ m - fi(r m , -)||o, J • First ' we use (lL22|) . (TOl) . (l2~9l) . (l2~6l) . (l2~5|) 
and (12.81). to obtain 



E 



r> VO J d 



< 



< 



x {Q,r m )( T ) [£m(T-;x,y) - *(t to - T\x,y)] W{T,y)dydi 
[fC m (r; x, y) - *(r m - r; x, y)] 2 dyefo^ dT 
51/ ( / / [ G ^."'- i + 1 ( x ^y) ~ *( Tm ~ T ',x,y)] 2 dydxj dr. 



dx 



Now, we introduce the splitting 
(4.23) 



/Km 
IJ 2 1 



where 



B i n '-=^2 (/ / [GA*, m -i+i{x,y) -^iTrn-n^x^)] 2 dydx) dr, 

Ja ( \J dJ d / 

B 2 f ill [^{T-,n-T^ 1 ;x,y)-^(T m -T;x,y)] 2 dydx\ dr. 
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By the definition of the Hilbert-Schmidt norm, we have 

rn og r / r p 

B^<AtJ2J2 ( / [ G A*, m -e+i(x,y)tpk(y) dy- I *(r m - Tt- X \x,y)<p k {y) dy 

og / m \ 

< E E Ar HA" 1 -**** - 5(r m - r,_x)^||^ 

fe=i \ #=i / 

CO / 1TL \ 

<E E Ar H Am "' + V' fc S(r m -n^UWln 
k=i \ e=i / 

oc / rn \ 

^E A ' E Ar||A^ fe -5(r,) efc ||L ■ 
fc=i \^=i / 

Let G [0, |). Using the deterministic error estimate f|4.5[) and (|2.10p . we obtain 

V^f< C(Ar) e ( f>||| £fc ||| 4e _ 2 j 



, fc=i 



(4.24) 



\ fe=i A »= 



Using, again, the definition of the Hilbert-Schmidt norm we have 

c>o rn 

B ™ = E E / - T ^ k - S ^ - T)<PX,n dr 

k=l £=1 ^ 
oo m p 

= E A fcE / \\S{T m - n-i)e k - S{T m -r)e k \\l D dT 

k=l t=l 

Observing that S{t)e k = e ~ A »( A *~")* e k for t > 0, (|4T25) yields 



(4.25) 



co m 



6 2 n = E A ' E 



(4.26) 



where 



fc=i £=i" /i « 
oo m 

= E A " E 

fc=i £=i" /i « 



k m 

^ 2 m i -E A ' E 

B. 



A» W D 



-(At-/*Aj)(r m -r < _ 1 ) _ e -(A£-M A=)(r m -T) 



1 _ g-l^-Mt) 



1 - e 



-M-/iA2)(T-T<-l) 



dr. 



oo m 



J 2,2 



E A " E 

k = K+l i = \ 



1 _ e -(At-^A=)(T-T^_!) 



dr. 
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First, we estimate B™ x and as follows 

oo 

S™2< E A2(l-e- A ^ A "-^ A ^) 2 



(4.27) 



OO 

E 1 

fc=re+l 



,-2(At-MA 2 fc )(r m -r) dT 



2,,2 



— 2(1+2k) 



c A fc ' 



fc=l 



(4.28) 



fe=l £=1 

k rn „ 

j i n i *■* A (i 



(It 



fc=i f=i ' 

\2 



< C (At) 

Finally, we combine (j4~27| . (|4T2"gj) and ([XT7]) . to obtain 



(4.29) V S 2 l < C (l + (Ar)* + (Ar)3j (At)*. 

The estimate (j4~2Tj) follows by (|433j) . ([4T24|) and (j4~29|) . □ 

5. Convergence of the Fully-Discrete Approximations 



□ 



To get an error estimate for the fully-discrete approximations of u defined by (|1.8I) - (|1.9[) . we proceed 
by comparing them with their time-discrete approximations defined by (|4.ip - (|4.2p and using a discrete 
Duhamel principle (cf. [T4 ] , [2"2" ] V 

5.1. The Deterministic Case. The Backward Euler finite element approximations of the solution to 
(|1.5|) are defined as follows: first, set 

(5.1) Wl ~ P h w , 

and then, for m = 1, . . . , M, find W r h n £ Ml such that 



WV 



ArA fl , h Wr = > 



(5.2) W? 

which is possible when ji 2 At < 4. 

Next, we derive a discrete in time L 2 (L 2 ) estimate for the error approximating the Backward Euler 
time-discrete approximations of the solution to (jl.5|) defined in (|4.3|) - (|4.4p . by the Backward Euler finite 
element approximations defined in (|5.1I) - (I5.2[) . The main difference with the case ^ = which has been 
considered in [14] . is that, our assumption (11.21) on /i, can not ensure the coerciveness of the discrete 
elliptic operator A B h . 

Theorem 5.1. Let r = 2 or 3, w be the solution to the problem (|1.5p . (VF m )^ =0 be the time-discrete 
approximations ofw defined in (|4.3[) - (|4.4[) . and (W / ™)*j =0 C M 7 h be the fully-discrete approximations of w 
defined in (|5.1[) - (|5.2p . Also, we assume that /j 2 At < j. If wq G H 2 (_D), then, there exists a nonnegative 
constant ci, independent of h and At , such that 



(5.3) 



E Ar \\ wr ' 



w, 



h \\0,D 



< ca h^ e WwoU^.eu V0€ [0,1], 
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where 
(5.4) 



L(r) r- 



2 if r = 2 
4 if r = 3 



and £*(r,0) := (r + 1)0-2. 



Proof. The error estimate (|5.3j) follows by interpolation, after showing that holds for 6 — and = 1. 
In the sequel, we will use the symbol C to denote a generic constant that is independent of At and h, 
and may changes value from one line to the other. 

Let E m := WJ? ~ W rn for m = 0, . . . , M. First, use (JO]) and gU) to obtain 



w? - wr 



At A B h W™ = At fj 2 W } 



2 TI/m 



(5.5) 

(5.6) W m -W rn ~ 1 +ATA B W m = ATfi 2 W m 

for m = 1, . . . , M. Then, combine (|5.5|) and (|5.6[) . to get the following error equation 

(5.7) f B , h {E m - E m ^) + At E m = At f f B . h E m - At (f B - T B ,„)A B W"\ to = 1, . 
Taking the L 2 (£>) —inner product with E m of both sides of (15.71) . it follows that 

7 B .„(£ m - E m -\E m ) + A T \\E m \\ 2 D = At /j 2 j B h (E m 7 E rn ) 

- At ((f B - T B ,„)A B W m , £™) 0>D , m = 1, . 
from which, after using (|2.1ip , we conclude that 

%, h (E m ,E m ) + At \\E m \\lo < l B , h {E m -\E m - x ) + 2 At M 2 %, h (E m , E m ) 

+ AT\\(f B -f Bth )A B W m \\l D , m = 1,...,M. 
Since 2At^ 2 < 1, (JEHJ yields 



,,M, 



(5.8) 



(5.9) 



1b.u(E 



m— 1 



AT||(T B -r B . h )A B W 



m||2 

0,D 



for to = 1, . . . , Af . Applying a simple induction argument based on (|5.8p and then using that 4At^ 2 < 1, 
we get 



(5.10) 



max %, h (E m ,E m ) < C 

< m < M 



%, h (E ,E ) + ATj2\\(Ts~T B , h )A B WX,n 

e=i 

Summing with respect to to from 1 up to M, using dHTTUJ) and observing that T BJl E° = 0, 

M A-/ 

(5-11) J2 AT W Em Wo-o< C J2 A T\\(T B -T B , h )A B W m \\t. D - 

m—1 m— 1 

Let r = 3. Then, by (12.22[) , (|5.1ip and the Poincare-Friedrich inequality, we obtain 

i i 

/ M \ 2 / M \ 2 

£ArpHlL> <Cfc 4 J] AT||A B iy m || 2 J 



gives 



(5.12) 



<C/i 4 



£ At ( ||9^ m |lo. D + ll^Vm 2 , + II^WHI„ 2 . . 



Taking the L 2 (D)~ inner product of (|4.4p with d 4 W m and then integrating by parts, we obtain 
(5.13) (d 2 W rn -d 2 W m - 1 ,d 2 W m ) 0tD + At ||d 4 W m ||£ D + Ai At (<9 2 tU m , d A W rn ) a . D = 0, to = 1, . . . ,M. 
Using (|2.11[) . (|5.13l) and the Cauchy-Schwarz inequality we obtain 

\\d 2 W m \\l B + 2 At ||d 4 WHIo, a < l|9 2 ^ m - 1 || 2 B + 2 M At ^W^X^ ||d 4 I¥ m || , D , to = 1, . . . ,M, 
which, after using the geometric mean inequality, yields 
(5.14) 



\d 2 W m \\Z D + Ar\\d*W m \\l D < \\d 2 W 



2W"-1||2 
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2 i|Q2 T i/m||2 



At fi \\d 2 W 



m = 1, 



.AT. 



Since 2^ 2 At < 1, from (|5.14j) follows that 

< e 2 M 2 Ar|| a 2 w m-l||2 Dj jjl — 1, . . . , M, 

from which, applying a simple induction argument, we conclude that 

(5.15) max \\d 2 W m \\ 2 , D < C |K|| 2 , D - 

0<m<M 

Next, sum both side of f|5.14D with respect to m, from 1 up to M, and use (|5.15p to conclude that 

M 

(5.16) Y, ^\\d±W m \\l D < C\\w Q \\l D . 

m— 1 

Taking the L 2 (D)~ inner product of (|4.4[) with d 2 W m , and then integrating by parts, it follows that 

(5.17) (dW m - dW m -\ dW m ) , D + At \\d 3 W m \\l D +ptAr (dW m , d 3 W rn ) a , D = 0, m = 1, . . . , M. 
Using (|2.11[) . (|5.17l) . the Cauchy-Schwarz inequality and the geometric mean inequality, we obtain 

\\dW m \\l D + At \\d 3 W m \\l D < WdW^fo, o + At fi 2 \\dW m \\l D , m = 1, . . . , M. 
Since 2 /j, 2 At < 1, proceeding as in obtaining (|5.15j) and ()5 . 1 6|) from (|5.14|) . we arrive at 



(5.18) max \\dW m \\l D + E Ar P^Wlo < C KIlL, 

0<m<M z — * 
m— 1 

Thus, combining ([S"T2]l . (pUB) . ([BTTB]) . (f5TT8|) and ([2l?]) . we obtain 

(5.19) (|f>r||£ m |lo, D J < Cft 4 || W o||^. 

Let r = 2. Then, by (|2.22[) , (|5.11|> and the Poincare-Friedrich inequality, we obtain 

i i 

(M \ 2 / Af \ 2 

<c/i 2 E Ar (ii^ m iio, D + ll^ m n 2 D ; 

_ m— 1 

Combining, now, (|5.20p . (|5.18p and (|2.3p . we obtain 

(5.21) (^jr Ar||£™|| 2 , D ) <Ch 2 \\wo\W. 

Thus, relations (|5~T9|) and flOT|) yield flO) and flO} for = 1. 
Since ^i 2 Ar < 1, using (|5.5|l . we have 

T s ,h(Wr " + At W h = A ^ T B , h W^, m = 1, . . . , M, 

from which, after taking the L 2 (D)— inner product with W™, we obtain 

(5.22) 7b,„ (WT* " ^r 1 , WT)o.x> + Ar ||^||J D = Ar ^ ^(WT. WT), m = 1, . . . ,M. 
Then we combine (|5.22p with p. lip to have 

(5.23) {l-2AT^ B , h {W^,W?) + 2AT\\wnlo< WW^" 1 , W^" 1 ), m = l,...,M. 
Since 4^ 2 Ar < 1, (|S~23l) yields that 

< e VAT Tb^OC" 1 , V^™" 1 ), m = 1, . . . , M, 
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from which, applying a simple induction argument, we conclude that 
(5-24) max %, h (WfT, WJT) < Cj B , h (W%, <). 

< ra < m 

Summing with respect to m from 1 up to M, and using (|5.24p , (15. 23ft gives 

M 

(5-25) At J2 WWnln < C%, h (W%, W%) , D . 

Finally, using f^5|) . P^gj) and (pEI|) we obtain 

^ Ar||Wr||= B <C(f B , h «;o,^o)o, D 



(5.26) m_1 



<C|ko||! 

<ciW2, 



2,D 

2 



Finally, combine (15T2T))) with (j4~T^]) to get 

^Ar||W ro -Wrilo, D J < C\\woU-,, 

which yields Q and (JEl for 6> = 0. □ 

5.2. The Stochastic Case. Our first step is to show the existence of a Green function for the solution 
operator of a discrete elliptic problem. 

Lemma 5.1. Let r = 2 or 3, e > with /i 2 e < 4, / £ L 2 (D) and iph G M£ suc/i t/iaf 

(5.27) */j h + eA B:h Tp h = P h f. 

Then there exists a function A e fi £ H 2 (D x D) such that A t) h \d(DxD) — an d 

(5.28) M*)= I A h , e (x,y)f(y)dy Viefl 



and Ah,e(x, y) = Ah,e(y, x) for x,y £ D. 

Proof. Let S e> h ■ x M 7 h — > R be the inner product on M£ given by 

&t,h{<t>, X) ■= e (A B ,h0, X)o,d + (0, x)o,n 



We can construct a basis (xj)"=i of M£ which is L 2 (D)— orthonormal, i.e., (Xi,Xj)o,n = 5ij for i,j = 
and c^-orthogonal, i.e., there exist {K,h,t)7=i c (0, +oo) such that 5e,h{XhXj) = K,h,i5%j 
for i,j = 1, . . . , n/, (see Section 8.7 in [9]). Thus, there are (/fj)J=i C R such that t/j/j = J^i=i Mj Xji an d 
(|5.27|) is equivalent to //, = A - (/, x<)°,d for i = 1, . . . ,71;,. Finally, we obtain (|5.28[) with Ah. t {x.y) = 

Our second step is to compare, in a discrete in time L^°(Lp(L 2 )) norm, the Backward Euler time- 
discrete approximations of u with the Backward Euler finite element approximations of u. 

Proposition 5.2. Let r = 2 or 3, u be the solution of the problem (|1.6[) . (C^™)m=o ^ e ^ e Backward 
Euler finite element approximations of u defined in (|1.8[) - f)1.9[) . and (C^ m )m=o i/ie Backward Euler 
time-discrete approximations of u defined in (|4.1j) - ()4.2|) . ^4Zso, we assume that \j? At < v. Then, there 
exists a nonnegative constant C2, independent of Ax, At, h and At, such that 

(5.29) max (e \\\UJ? - U m f Q JlY < c 2 e~^ h v ^~\ Ve G (0,i/(r)], 

0<m<Af \ L" iiu.jjj / 



where 
(5.30) 



if r = 2 
if r = 3 



Proof. Let I : L 2 (D) — > L 2 (D) be the identity operator and : L 2 (D) —> M" h be the inverse discrete 
elliptic operator given by A^ := (I + At A B>h )~ Ph, having a Green function G A — Ah t Ar according to 
Lemma r5.il and taking into account that /i 2 At < 4. Also, wc define an operator $h ■ L 2 (D) — > MJ^ by 
($ft/)(x) := J n G* h (x,y) f{y) dy for / e L 2 {D) and x e D, where G^ h (x,y) = -d y G Ah {x,y). Then, we 
have that Ahf — §hf f° r au / € H 1 (D). Also, for I £ N, we denote by G A t < the Green function of 
A e h &h- In the sequel, we will use the symbol C to denote a generic constant that is independent of At, 
Ax, h and At, and may changes value from one line to the other. 
Applying, an induction argument, from (I1.9P we conclude that 



Uh = E / K~ j ®hW(T, •) dT, m = l, 

3 = 1 JA 3 



,M, 



which is written, equivalently, as follows: 



(5.31) 



J D 



V h . m {T;x,y)W{T,y)dydT VxED, m = l,...,M, 



where V h>m (r;x,y) := £™ 1 X a .(t) G A (a;, y) Vr G [0,T], Vx,y e £>. Using flL22]), ([OT]), the 

Ito-isometry property of the stochastic integral, (|2.5[) and the Cauchy-Schwarz inequality, we get 



E 



\\u m -u?\\t D < 



< 



[k. m (T;x,y) - V h>m {T;x,y)\ dydxjdT 
J2 [ \\A m -^-A"r 3 <£ h \\l s dT, m=l,...,M, 



where A and $ are the operators defined in the proof of Theorem 14.21 Now, we use the definition of the 
Hilbert-Schmidt norm and the deterministic error estimate (|5.3[) . to obtain 



E 



\u m ~uni»] < E Ar 



i=i 

OO 

^ E 



< 



k=l L 1=1 

oo 

fc=l 



k=l 

EAr||AV fe -Ai^|| 2 D 
E Ar||A^ fc -Af l£fc || 2 



fk\\L 



<Ch 2i ^ e J2\t\\e k \\l Mr , e) , m = l,...,M, V0e[O,l]. 



fc=i 



Thus, we arrive at 



(5.32) max 

Km<M 



(e[|| 



m II 2 
ft, llo.D 



< c//*« e E A 



, V0e[O,l] 



. k=l 



It is easily seen that the series in the right hand side of (|5.32[) convergences iff v(r) > i*(r)d. Thus, 
setting e = v(r) — l+(r)6, requiring e <= (0, v{r)\, and combining (|5.32p and (|2.10[) . we arrive at the 
estimate (j!T2T)]) . □ 

The available error estimates allow us to conclude a discrete in time Lf'(L 2 p (L 2 )) convergence of the 
Backward Euler fully-discrete approximations of u. 
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Theorem 5.3. Let r = 2 or 3, i/(r) be defined by (|5.30|) . u be the solution of problem (|1.6[) . and (UVf)^_ Q 
be the Backward Euler finite element approximations of u constructed by (ll.8[) - (|1.9[) . Then, there exists 
a nonnegative constant C, independent of h, At, At and Ax, such that: if /i 2 At < j, then 



max 

0<m<M 



{E[\\Ujr-u(r m ,-)\\l D ]} 



< C 



w*(Ar,ei) Ars~ £l + e 2 



forall ei € (0, §] and e 2 G (0,i/(r)], w/iere w»(Ar,ei) := e x 2 + (Ar) ei (l + (At) a + (Ar)i)i 
Proof. The estimate is a simple consequence of the error bounds (|5.29l) and (|4.21[) . 



□ 
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Appendix A. 

Let t > and /i/- := A|(A| — //) for k e N. First, we recal that S(t)wo = Y^k^i e_Mfc ' i w O, s k)o,n £k for 
t > 0, and set S(i)wo = * S(t)wo for t > 0. Next, follow Chapter 3 in [3T], to obtain 



\d e t S(t)w \\l P = XX (dfS(t)w ,s k )l o 

k=l 

= f>f (M fe +M 2 r (S(t^o,e k f OD 

k=l 

oo 



2 

0,-D' 



fc=l 



which yields 

oo 

(A.l) \\d e t S(t)w \\l P < £ Af +u h-^ (w ,e k )l D , 

k=l 

~ / 2 \ 2f 

where := ( 1 + ^ + ^ J ■ Now, use (|A.1|) . to have 

(r-O^l^rKlLA-^G^ £ A f +4 ^ 2 « ( / [\i(r ~ t a f dr) (w ,e k )l D 

k=l Jt * 



fe=i 







. />oo x 00 

XJo 1 k=i 

which yields 

(A.2) [\r- t a f \\dfS(r)w \\l P dr < C 0Afl \\w \\l P+ u- 2fl - 2 , 



where Cp Alt = / °° e~ x dx. Observing that dfS(t)w = 4 E„=o O M 2( '" m) and 
using (|A.2|) . we conclude that 

/ (r-taf \\d e t S(r)wo\\l P dT < e 2 ^ T Cm,* Y, IkoIlL^™-^ 

J t„. n 



m=0 



which yields HHg) with C^.^t = Cp. £ ^ e 2 ^ T £. □ 
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